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Abstract. The nature and the location of the lenses discovered in the microlensing surveys done so far towards the 
LMC remain unclear. Motivated by these questions we compute the optical depth and particularly the number 
of expected events for self-lensing for both the MACHO and ER0S2 observations. We calculate these quantities 
also for other possible lens populations such as thin and thick disk and galactic spheroid. Moreover, we estimate 
for each of these components the corresponding average event duration and mean mass using the mass moment 
method. By comparing the theoretical quantities with the values of the observed events it is possible to put some 
constraints on the location and the nature of the MACHOs. Clearly, given the large uncertainties and the few 
events at disposal it is not possible to draw sharp conclusions, nevertheless we find that certainly at least 3-4 
MACHO events are due to lenses in LMC, which are most probably low mass stars, but that hardly all events 
can be due to self-lensing. This conclusions is even stronger when considering the EROS2 events, due to their 
spatial distribution. The most plausible solution is that the events observed so far are due to lenses belonging to 
different intervening populations: low mass stars in the LMC, in the thick disk, in the spheroid and possibly some 
true MACHOs in the halo. 
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1. Introduction 

The location and the nature of the microlensing events 
found so far towards the Large Magellanic Cloud (LMC) 
is still a matter of controversy. The MACHO collabora- 
tion found 13 to 17 events in 5.7 years of observations, 
with a mass for the lenses estimated to be in the range 
0.15 — 0.9 Mq assuming a standard spherical Galactic halo 
(Alcook et al. 2000a) and derived an optical depth of 
T = 1.21q3 X 10~^. An analysis of the spatial distribution 
of the events as well as a detailed study of the source star 
location done on HST images (Alcook et al. 2001a) are 
both consistent with an extended lens distribution such 
as the Milky Way halo, however an LMC distribution is 
only slightly disfavoured. Thus this test is not conclusive 
given the few events at disposal. The ER0S2 collabora- 
tion (Milsztajn & Lasserre 2001) announced the discovery 
of 4 events, with an estimated average mass of 0.2 Mq, 
based on three years of observation but monitoring about 
twice as much stars as the MACHO collaboration. The 
MACHO collaboration monitored primarily 15 deg^ in the 
central part of the LMC, whereas the EROS experiment 
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covers a larger solid angle of 64 deg^ but in less crowded 
fields. The EROS microlensing rate should thus be less 
affected by self-lensing. This might be the reason for the 
fewer events seen by EROS as compared to the MACHO 
experiment. 

It has been argued that some if not all the events could 
be due to LMC self-lensing. Indeed, several authors have 
already studied in detail self-lensing by focusing, however, 
mainly on the value of the microlensing optical depth. 
Sahu (1994) and Wu (1994) suggested that self-lensing 
within the LMC could explain the observed optical depth. 
This claim has been questioned by further studies (Gould 
1995; Alcock et al. 1997a). A major problem is the un- 
certainty related to a precise knowledge of the shape and 
the total mass of the LMC. Indeed, for instance a dis- 
aligned bar, as suggested by Weinberg, could increase 
the self-lensing optical depth (Weinberg 2000; Evans & 
Kerins 2000; Zhao & Evans 2000) as well as assuming the 
LMC being much more extended along the line of sight 
(Aubourg et al. 1999). These authors find then an optical 
depth for self-lensing of 0.5 — 1.5 x 10^^, thus comparable 
with the measured value by the MACHO collaboration. 
Other authors, instead, find lower values for the optical 
depth in the range 0.5 — 8.0 x 10~^ (see for instance Gyuk 
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et al. 2000). The discrepancy are of course due to the dif- 
ferent adopted models for the shape of the LMC, which 
arc based on various arguments and by weighing differ- 
ently the observation on the LMC star distribution which 
are available today. 

Wc notice that a direct comparison of the theoretical 
values for the optical depth to self-lensing and the mea- 
sured one is not straightforward. Indeed, the theoretical 
values arc computed for a given line of sight and then 
to compare with the measured value one takes the aver- 
age over the lines of sight corresponding to the various 
monitored fields (Gyuk ct al. 2000). On the other hand 
the measured optical depth is computed assuming that 
it stays constant over the monitored fields of the LMC. 
Such a procedure is certainly adequate if the lenses are in 
the halo and thus their optical depth does practically not 
vary on the size of the LMC, however, this is no longer 
true when dealing with self-lensing. 

Some of the events found by the MACHO team are 
most probably due to self-lensing: the event MACHO- 
LMC-9 is a double lens with caustic crossing (Alcock ct 
al. 2000b) and its proper motion is very low, thus favour- 
ing an interpretation as a double lens within the LMC; 
the source star for the event MACHO-LMC-14 is double 
(Alcock et al. 2001b) and this has allowed to conclude that 
the lens is most probably in the LMC. The expected LMC 
self-lensing optical depth due to these two events has been 
estimated to lie within the range 1.1 — 1.8x 10~* (Alcock et 
al. 2001b), which is still below the expected optical depth 
for self-lensing even when considering models giving low 
values for the optical depth. 

The event LMC-5 is due to a disk lens (Alcock et al. 
2001c) and indeed the lens has even been observed with 
the HST. The lens mass is either ~ 0.04 Mq or in the 
range 0.095 — 0.13 Mq, so that it is a true brown dwarf 
or a M4-5V spectral type low mass star. The other stars 
which have been microlensed were also observed but no 
lens could be detected, thus implying that the lens cannot 
be a disk star but has to be cither a true halo object or a 
faint star or brown dwarf in the LMC itself. 

Some work has also been done on the Small Magellanic 
Clouds (SMC) where only two microlensing events have 
been found up to now (Alcock et al. 1997b, 1999; 
Palanque-Delabrouille et al. 1998). One was a resolved 
binary event, allowing the determination of the lens dis- 
tance (Alcock ct al. 1999; Afonso ct al. 1998; Albrow et al. 
1999), which most probably resides in the SMC itself and 
thus clearly being due to self-lensing. The other event is of 
long duration and a detailed analysis seems also to favour 
a self-lensing interpretation (Palanque-Delabrouille et al. 
1998). It has been argued that the SMC self-lensing opti- 
cal depth should be higher than the corresponding LMC 
value since the SMC is tidally disrupted, due to its in- 
teraction with the Milky Way and the LMC, and is thus 
quite elongated along the line of sight to the Milky Way 
(Caldwell & Coulson 1986; Welch et al. 1987). 

Thus up to now the question of the location of the 
observed MACHO events is unsolved and still subject to 



discussion. Clearly, with much more events at disposal one 
might solve this problem by looking for instance at their 
spatial distribution. However, since the MACHO collab- 
oration data taking stopped at the end of 1999 and the 
EROS experiment is still underway, but will hardly lead 
to a large amount of events in the next few years, there is 
not much hope to have substantial more data at disposal 
within the next few years and it is thus of importance 
to explore further ways which can give insight into the 
MACHO location using only the already existing data. 
This is the main aim of this paper. 

As a first point we calculate the optical depth, the av- 
erage duration and particularly the number of expected 
events for self-lensing. For the latter quantity we take for 
the number of monitored stars and for the exposure time 
the values corresponding to the MACHO and EROS ex- 
periments, respectively. Moreover, we compute the same 
quantities for MACHOs located in a thin or thick disk or 
in a spheroid around our Galaxy. As a main point we com- 
pute following the mass moment method (De Rujula et al. 
1991, 1992) the average mass of the lenses for the observed 
events under the assumption that the lenses are located 
in the LMC, the halo or in one of the Galactic compo- 
nents. This allows to see whether the assumption of hav- 
ing the lenses in the LMC itself or in one of the Galactic 
components leads to consistent values for the masses. Our 
analysis indicates that most probably not all lenses orig- 
inate from the same population but are due to different 
ones. Therefore, some of the conclusions reported so far in 
the literature, especially on the average mass and thus on 
the nature of the lenses, have to be taken with caution. 
Finally, we give an estimate of the fraction of the local 
dark mass density detected in form of MACHOs in the 
Galactic halo. 

The paper is organized as follows: in Sect. 2 wc give 
the density profiles of the stellar distributions for the dif- 
ferent components of the LMC, which will then be used 
in the following. In Sect. 3 wc present the different equa- 
tions to compute the optical depth, the microlensing rate 
and the mass moments we shall use. In Sect. 4 wc report 
the results for the LMC self-lensing, whereas in Sect. 5 
we present the values we find for the various Galactic lens 
populations: thin and thick disks, spheroid and halo. The 
mass of the lenses is discussed in Sect. 6. In Sect. 7, by 
using different models, we determine the fraction of the lo- 
cal dark mass density detected in form of MACHOs in the 
Galactic halo. In Sect. 8 we analyse the asymmetry in the 
spatial distribution of the observed events. We conclude 
in Sect. 9 with a summary of our results. 

2. LMC morphology, mass and density profiles 

In the last years the puzzle of the real nature of the LMC is 
showing up progressively with respect to the first picture 
of de Vaucouleurs (1957). The unexpected scale length 
substantially larger than most irregular galaxies led to re- 
vise the LMC shape (Bothun & Thompson 1988). Kim et 
al. (1998), emphasizing the spiral structure of the HI gas of 
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LMC, confirmed the opinion that LMC has to be consid- 
ered, from a morphological and dynamical point of view, 
more similar to a dwarf spiral than an irregular galaxy. 
However, LMC does not present a well-defined center, 
and there are still uncertainties and discordant measures 
on many important parameters. For this reason we will 
use different models and vary some parameters in specific 
ranges such as to get reasonable ranges for our theoretical 
results. 

The LMC conventional features are a disk and an off- 
centered bar lying in the same plane. The significant offset 
of the bar in respect of the disk (about 1 kpc) suggests a 
LMC dominated by dark matter (Cioni et al. 2000), in 
agreement with the hypothesis of Still (1999) that LMC 
is a fast rotator dwarf galaxy. There are many lines of 
evidence supporting the planar geometry of LMC (van 
der Marel & Cioni 2001). Nevertheless, several authors 
suggested that some LMC populations can not reside in 
the disk plane (Luks & Rohlfs 1992; Zaritsky & Lin 1997; 
Zaritsky et al. 1999; Weinberg & Nikolaev 2001) while 
Zhao and Evans (2000) proposed a misaligned offset bar 
model. 

The surface-brightness profile of the LMC is strongly 
exponential along the bar and a truncation becomes ap- 
parent at a scale length of 1.65 kpc (Bothun & Thompson 
1988). The disk seems circular, thin, flat and seen nearly 
face-on with the east side closer to us than the west side. 
Its plane is tilted to the plane of the sky by an angle of 
roughly 30° — 45° with a position angle of (j>d ~ 170° 
(Westerlund 1997). Following Gyuk et al. (2000), we 
placed the disk center at the kinematic center of the HI 
gas J2000 (Zd,o, bd.o) = (279.7°, -33.5°) (Kim et al. 1998). 

Thanks to the recent DENIS and 2MASS surveys, 
there are strong indications that the shape of the LMC 
disk is not circular, but elliptical with a small vertical 
scale height (< 0.5 kpc) (van der Marel & Cioni 2001). 
The study of near-IR surveys yields an intrinsic elliptic- 
ity of 0.31 in the outer parts of the LMC (van der Marel 
2001). This elongation, considerably larger than typical 
for disk galaxies, is in the direction of the Galactic center 
and perpendicular to the Magellanic stream, and is likely 
caused by the tidal force of our Galaxy. Very recently, van 
der Marel et al. (2002) pointed out that the distribution 
of neutral gas is not a good tracer, and thus leads to an 
incorrect LMC model. Instead, using the carbon star data, 
they provided an accurate measurement of the dynamical 
center of the stars, which turns out to be consistent with 
the center of the bar and with the position of the center 
of the outer isophotes of the LMC. Thus, the disk has a 
center that coincides with the center of the bar, has an 
elliptical shape and is thick and fiared (van der Marel et 
al. 2002). We will discuss a new LMC model that takes 
completely into account this revised picture of the struc- 
ture and dynamics of the LMC in a forthcoming paper. 
However, we do not expect the results of this paper to 
be substantially modified and thus also our main conclu- 
sions, since the reported total mass of the LMC is well in 
accordance with the values we adopted. 



The bar has a size of roughly 3° x 1°, a position 
angle of 0b — 120° and its center is located at J2000 
(^b,o,6b,o) = (280.5°, -32.8°) (NEDQ data). The bar has 
boxy contours and sharp edges, even if new morphologi- 
cal studies, based on colour-magnitude diagrams extracted 
from DENIS catalogue, show a bar extending over about 
4° as well as an elongation at its borders that suggests the 
presence of spiral arms (Cioni et al. 2000). 

There is a great debate also about the LMC mass. This 
topic is discussed in detail in Gyuk et al. (2000). Here, 
according to several authors (e.g. Sahu 1994; Gyuk et al. 
2000), we adopt A/bar + Mdisk in the range 3 - 6 x IO^Mq 

with 1 < Mdisk/Mbar < 5. 

Another important parameter is the velocity disper- 
sion: for the LMC stars, we adopted a velocity dispersion 
in the range of 15—35 km/sec, as it seems that there are no 
LMC stellar populations with a velocity dispersion greater 
than 33 km/sec (Hughes et al. 1991). 

In the self-lensing framework, we must consider differ- 
ent source/lens geometries, e.g. sources and lenses in the 
disk or in the bar, otherwise sources in the disk, lenses 
in the bar and viceversa. Due to the considerable uncer- 
tainties of the various LMC components, we will consider 
two different models for the description of the LMC disk 
and bar populations, model 1 and model 2. Each of the 
models are then further subdivided into two classes a and 
b according to the set of the adopted parameters. 

• Model 1. Following Gyuk et al. (2000), we consider 
-Mbar = 1/4 Afdisk, and adopt a circular LMC star den- 
sity profile to describe the spatial distribution for the 
disk population, which will be used for both the lenses 
and the sources. We can write their density as a double 
exponential profile 



Pd 



47rzd-R3 



(1) 



where i?d = 1.6 kpc is the radial scale length, Zd = 
0.3 kpc the vertical scale height and = 2.4 — 4.8 x 
10^ Mq the disk mass. R and z are the cylindrical 
coordinates of the lens/source position. For the disk 
inclination i we take a value in the range 30° to 40°. 
For the bar lens and source populations we use a tri- 
axial gaussian density profile to describe their distri- 
bution as follows 



Pb 



Mb 



{2n f^^ XhVhZb 



(2) 



where Mb — 0.6 — 1.2 x 10^ Mq is the bar mass, x, y 
and z are the coordinates along the principal axes of 
the bar. Due to the uncertainties we used two different 
sets for the scale lengths along the axes: 

— Model la with Xb = 1.0 kpc and yh ~ ~ 0.3 kpc; 

— Model 16 with xi, — 1.51 kpc and yb = ^b = 
0.44 kpc. 



^ Nasa/ipac Extragalactic Database 
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All other parameters in models la and lb are equal 
and have values as mentioned above. 

Model 2. Following the results on LMC morphology 
by van der Marel (2001), we consider an elliptical disk 
with an inclination of 35° and a line-of-nodes position 
angle equal to 122.5°. Zhao and Evans (2000) used a 
LMC lensable mass in the range 2.5— 5.5x IO^Mq, with 
the bar making up between 25% — 50% of it, suggesting 
thus that bar and disk could even have the same mass. 
This is supported by the fact that half of the LMC total 
number of stars is in the bar and this factor increases 
for younger objects (Cioni et al. 2000). We take a bar 
size of 4° X 1° (Cioni et al. 2000). 
Taking into account all these facts, we model the disk 
as 



Pd 



AnxdydZd 



(3) 



where x, y and z are the coordinates along the prin- 
cipal axes of the disk, and the scale lengths are — 
1.6kpc, yd = 0.7 Xd ~ 1.12kpc and Zd — 0.3 kpc. 
For the bar density we use a quartic exponential pro- 
file, to take into account its boxy shape: 



Pb 



^3/2 r (I) 



(4) 



where the scale lengths are Xh = 2.1 kpc, i?b = 
0.58 kpc, r (I) is the Euler gamma function and R 
is the radial coordinate in the orthogonal plane with 
respect to the bar principal axis x. In order to take 
into account the uncertainties in the LMC bar and 
disk masses, we again considered two cases for model 
2: 

- Model 2a with Mbar = l/4Mdisk and A/bar + 
Afdisk = 3-6x 109 Mq; 

- Model 2b with Afbar = Afdisk and Afbar + A/disk = 
3-5 X IO^Mq. 

All other parameters of both model 2a and 2b are equal 
and have values as mentioned above. 

With the help of the coordinate transformation given in 
Appendix A (see also Weinberg & Nikolaev 2000), we can 
easily rewrite R, x, y and z as a function of the variables 
observer-lens distance, Doi, observer-source distance, Dos, 
and of the parameters Dq, I, b, Iq, bo, (pd, (ph, i- Here 
Dq = 50 kpc (van der Marel et al. 2002) is the distance 
between the observer and the LMC, while I and b are the 
usual galactic coordinates. 

Again, following Gyuk et al. (2000) we will also con- 
sider a LMC halo contribution, for which we adopt a 
spherical model with density profile 



R^ 

a2 



(5) 



where a « 2 kpc is the LMC halo core radius (Gyuk et al. 
2000) and po the central density (see Table 0). We use an 



extreme model with a LMC halo mass of 6x 10^ Mq within 
8 kpc and a velocity dispersion tTh = 70/-\/2 km/sec, with 
a tidal truncation radius of 11 kpc (Gyuk et al. 2000). 

3. LMC self-lensing microlensing quantities 

3.1. Optical depth 

An important measurable quantity in a microlensing ex- 
periment is the optical depth r, which is defined to be the 
probability that at any time a random star is magnified by 
a lens by more than a factor of 1.34. For MACHOs in the 
Galactic halo we can in good approximation assume that 
all sources in the LMC are at the same distance Dos = Dq. 
The optical depth to gravitational microlensing is then de- 
fined as follows 



-Do 



Pi ■ 



Doi(Dos - Do 



Do 



■ dDoi, 



(6) 



where p\ denotes the mass density of the lenses, for in- 
stance MACHOs in the halo. 

The assumption that the sources lie all at the same dis- 
tance is no longer acceptable if we consider self-lensing. 
In this case we have to integrate not only on the dis- 
tance of the lenses but also on the distance of the 
sources. Moreover, the LMC extension can not be ne- 
glected, since the number density of the sources and lenses 
can change substantially with the distance (Kiraga & 
Paczynski 1994). Assuming that the number of detectable 
stars varies with the distance as Dof , where /3 is a param- 
eter that takes into account the magnitude limitations of 
the observation (a reasonable range for it is —3 < (3 < —1), 
the Eq. (g) becomes 



oo \ — 1 

Ps Dl+^^dD^, I X 



(7) 





Doi(Dos-Doi) ,n 
Pi aDoi 

^os 



where Ps denotes the mass density of the sources. 

Clearly for self-lensing the optical depth can vary sub- 
stantially with the position of the considered target field. 
Therefore, one has to compute the optical depth per unit 
surface and then perform an integration over the area 
which has been observed in order to compare the theo- 
retical value of the optical depth for self-lensing with the 
measured one reported in the literature. 

3.2. Microlensing event rate 

To estimate the number of microlensing events with a 
magnification above a certain threshold (where Ath — 
{ul^ + 2)/(wthVut\ -K4), Ah = 1-34 for wth = 1), we in- 
troduce the differential number of microlensing events 

dTVov = NJohsdT, (8) 

where N^, is the total number of monitored stars during the 
observation time tobs- dV is the differential rate at which 
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a single star is microlensed by a compact object and it is 
just equal to 

n{x)f (vi) d^x d^v 



dT = 



dt 



(9) 



where the numerator on the right hand side is the number 
of MACHOs with velocity in d^v = dvxdvydvz around vi, 
and located in a volume element d^x = dxdydz around 
the position x in the microlensing tube. n{x) is the lens 
number density and f{v\) is the velocity profile (VP) of 
the lenses, which is often assumed as having Gaussian 
shape characterized by a mean velocity and a velocity 
dispersion. Indeed, this is strictly true only for a singu- 
lar isothermal spherical model. However, the computation 
of the VPs (in particular along the line-of-sight) for non 
spherical shapes gets quite complicated as one has to solve 
the collisionless Boltzmann-Vlasov equation. Only for few 
specific models, so called power-law models (Evans & de 
Zeeuw 1994), an analytical solution has been found. In 
order to overcome this difficulty several authors (van der 
Marel & Franx 1993; Gerhard 1993) proposed to write the 
VPs along the line-of-sight as a new set of Gauss-Hermite 
moments, which arise from the expansion in terms of cor- 
responding orthogonal functions: the Gauss-Hermite series 



1 



V2 



N 



IT a- 



(10) 



J=3 



The zeroth-order term is a Gaussian and the higher or- 
der terms measure deviations from a Gaussian. The de- 
viations are quantified by the Gauss-Hermite moments 
/13, Hn'- the even coefficients /12/ quantify symmetric de- 
viations, while the odd coefficients /12/-1-1 anti-symmetric 
deviations. The Hj are the Hermite polynomials. Detailed 
studies show that the difference with respect to assum- 
ing a pure Gaussian shape, both for elliptical and spiral 
galaxies, is roughly 10 — 20% (van der Marel et al. 1994; 
Dehnen & Gerhard 1994). Moreover, the use of the Gauss- 
Hermite series introduces more extra parameters. On the 
other hand the uncertainties of the LMC parameters, as 
its total mass, are such that the theoretical estimates for 
microlensing quantities vary already widely for the allowed 
range of values. One has also to consider that the velocity 
dispersion of the lenses differs according to the consid- 
ered stellar population. This last fact is in any case diffi- 
cult to take into account when solving for a self-consistent 
solution in the framework of the collisionless Boltzmann- 
Vlasov equation. Thus for these reasons we prefer to re- 
strict to the lowest order of the VPs and, instead, vary 
the parameters, especially the velocity dispersion, in an 
enough large range. Therefore, following also Han & Gould 
(1995), for the VP we take a Maxwellian distribution with 
dispersion velocity cti, that is 



f{v,)d^v = 



1 



(11) 



It is convenient to write v\ = v\^~\-v^ and d'^v = d^VA_dv]i^ , 
since the component of the MACHO velocity parallel to 



the I.O.S., v\\ — VD^n does not enter in the description of 
the microlensing phenomenon. We will thus perform an 
integration of Eq. ( |ll| ) on this velocity component, such 
as to obtain the transversal velocity distribution (Jetzer 
1994): 



(12) 



Moreover, when writing the volume and the velocity ele- 
ments, we must also take into account the transverse ve- 
locities of the source star, i;sj_, and of the observer, VQ^^ 
(here we consider the observer co-moving with the Sun), 
that is the microlensing tube moves with a transverse ve- 
locity equal to 



, Dox\ , Col 



(13) 



So, we can write the volume element, d^x, as 

d^x = {vr± ■ n)dtdS = \vr±\ cos OdtdldDoi (14) 
= v,:± cos 0dtR^duthdadDo\, 

where Vr± — vi± — Vtj_ is the lens transverse velocity in 
the rest frame of the Galaxy, 



AGMq f DoiiD,,-D,i) 



(15) 



is the Einstein radius, while 9 is the angle between v^-^ 
and the normal, n, to the lateral superficial element, dS = 
dldDoi, of the microlensing tube, with dl = R-E,duthda 
being the cylindrical segment of the tube. Here and in the 
following, we use solar mass units defined as = M/Mq, 
where M is the lens mass. The velocity element, d'^v±, is 
given in cylindrical coordinates by 



d^v± — V,:±dVr±d9. 



(16) 



The number N of MACHOs inside a surface element of 
the microlensing tube becomes 

dN — n{x)f {vi±) v'^j_ cos dt i?E dut\i da dDd dvr±d9. {17) 

Therefore, the microlensing differential event rate is just 
(De Riijula et al. 1991, Griest 1991) 

rfP = — d/i/(t)i^)w^| cosO R^dutY^dadDoidv-c^dO. (18) 
d^, 



3.3. Mass distribution 

The determination of the mass distribution of the objects 
which act as lenses is one of the main aims of microlensing 
experiments. We assume, as usually, that the mass distri- 
bution of the lenses is independent of their position in the 
LMC or in the Galaxy (factorization hypothesis). So, the 
number density can be written as 

p)d,= ( P-) J^d, = ( ^) mD.,)d,, (19) 

dflj \ dfl J po.lens \ dfJ, J 
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where Hi{Do\) = pions/po.icns, Pions and po.ions are the den- 
sity and the local density of the lens population, respec- 
tively. The lens number density per unit mass, duQ/dfi, is 
normalized as follows: 



dno , Po.lens 



(20) 



where we will assume in the following for the lenses a 
minimal mass /imin = 0.1 or 0.01 and a maximal mass 
A'max = 10. If all lenses have the same mass /zq, we can 
describe the mass distribution by a delta function 



POJcn; 



-r- = 77 ~ Mo) 



(21) 



Otherwise, we can use a Salpeter-type initial mass func- 
tion (IMF) 



dno 
dfjL 



(22) 



where the classical value for a is 2.35 (Salpeter 1955). 
Even if this value overestimates the number of low-mass 
stars with M < 0.5 Mq, in general it describes very well 
stars with mass above ^ 0.3 — 0.5 Mq in different regions 
of the Milky Way and of the nearby Magellanic Clouds. 
Another choice, based upon counts of M dwarfs in the 
Galactic disk, is a equal to 0.56 for /i in the range 0.1 
to 0.59 and to 2.2 for n > 0.59 (Gould et al. 1997). In 
this case, the pre-factors Ci are fixed by the normalization 
condition 



Ci 



0.59 



,1-0.56 



dfl + C2 I fJ. 

0.59 



1-2.2 



dfl 



PO,lcns 



M, 







(23) 



We notice that, for these IMF slopes, /imin might even be 
lower than 0.1. Several alternative functional forms have 
been recently proposed for the low mass and brown dwarf 
stars range with a fimin — 0.01. Reid et al. (1999) indi- 
cated a value of a from 1 to 2, based on ultracool dwarfs 
discovered by DENIS and 2MASS surveys. By analyz- 
ing methane dwarfs, Herbst et al. (1999) suggested that 
a < 0.8 for disk brown dwarfs. Bejar et al. (2001), study- 
ing the substellar objects in the a Orionis young stellar 
cluster, found that the mass spectrum increases towards 
lower masses with an exponent a « 0.8 for 0.01 < ii < 0.2. 
Given all these uncertainties we will in the following per- 
form the various calculations using four different types of 
IMFs (see Sect. 4^), two with fimin = 0.1 and two with 

^min — 0.01. 

3.4. Mass and time moments 

Assuming that only the faint stars up to about /i^p — 1 
can contribute to microlensing events, the mass moment 
method (De Riijula et al. 1991; Grenacher et al. 1999) 
allows to extract information on the average mass of the 
lenses, relating the m-th mass moment 



£nitJ.)-j—tJ-"'dfi, 



(24) 



to the cumulative n-th time moment, constructed from the 
observations as follows 



< T" >= J2 ^" 

events 



(25) 



with m — {n + l)/2. T is the duration of a generic mi- 
crolensing event, that is 



Re 



Re 



\Vr±\ \Vij_-VtA 



(26) 



while Enifi) is the efficiency function (see De Riijula et 
al. 1991). In this way, with the knowledge of the event 
durations, we will be able to estimate the mean lens mass. 
The theoretical expression for the time moments is 



< T" >= 



£„(/i)r"d7Ve. 



(27) 



Clearly under the hypothesis (^ the Eq. (^ factorizes. 
The rate of microlensing events is given by: 




dVr 



2tv 





X / cos9d0 I duth / eol^j,)^^— — d^. 



Hi{Doi)dDoi / da (28) 



The 6 integration is between (— 7r/2,7r/2), since we are 
only interested in lenses entering the microlensing tube. 

As we already mentioned in the above section, since 
we are considering LMC self-lensing, the Eq. (^8|) has to 
be integrated not only over the distance of the lenses but 
also over the distance of the sources. 

Till now we have considered the most general situation. 
However, because the distance among lenses and sources 
in the self-lensing hypothesis is small when compared with 
the observer-lens or observer-source distance, we have that 
^ « 1 (Alves & Nelson 2000). Therefore, the transverse 



tube velocity, as defined in Eq. (13), is roughly equal to 
the source velocity, Vt± ~ Vs±, and the integral in da 
becomes trivial. Indeed, we have checked numerically that 
the results with or without this approximation differ very 
little, so that the induced error is negligible. 

Moreover, we must also take into account the fact that 
the source stars are not at rest but have a velocity Vs- 
Since we consider that both the lenses and sources are in 
the LMC, we can use another Maxwellian to describe the 
velocity distribution of the sources, with the same velocity 
dispersion (cs = ai = a) as for the lenses: 



F{vs±)dvs±d(p = -^e dvs±d(p, 



(29) 



where ip is the angle between and Vs±. Therefore, 
taking into account all these facts, the microlensing rate 
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becomes 



and the lens mean mass 



2ir 







X / Vs^_e dvs± / v^j^e \ ■^i ) dvrA_, (30) 

^0 "'0 



where 



Das 



Do 



■Hi{Doi)dD. 



ol 



X H,{Dos)Dlt^f'dDos, (31) 

with Ha = H,{Dos)Dl+^l'dDos and iIs(i?os) - 
Ps/pa.s- Then, adopting for the constant /3 the value 
P = —1 (which should be reasonable when considering 
main sequence stars), the general expression for the n-th 
time moment is: 



<T" > = y"e„(/i)T"diVev = 



\Vr±J 
<■ P- > 



(32) 



with 



jn 



27r 



ttHo 



OO pOO 



V c2 



/o {2mV2) r/ir,2-"e-2^''-''^'d?7id?72,(33) 

10 Jo 

where ryi — Vs_\_/ct, r/2 — Wix/t, /q is the modified Bessel 
function of order 0, fc = UminA^^iobs, and 



Das 



Do. 



Hi{Doi)dD. 



ol 



X Hs{Dos)dDo 



(34) 



We are able now to write the time moments < >, 
<T^ > and < T^^ >: 



< > 



< Ti > 



<f,2 >j(i 2), 



AGM, 



,0 



r^<M'>7(l), 



<T-i > =</>7(0) 
where 

7(l/2) = ^/."/i/^ 

7(1) = ^/.^/^, 



7(0) 



4A: 



1 rO 



(35) 

(36) 
(37) 

(38) 
(39) 
(40) 



The microlensing quantities of interest are the average 
event duration 



<T >-- 



< TO > 



< > _ ugmq < > 7 (1) 



c2 </,^>7(l)' 



(41) 



< M lii^ ^ Ali^. (42) 



< ^0 > y c2 



< T-1 > 7(1)' 



For the mean mass, according to Eq. ([25|), the time mo- 
ments < > and < > are calculated using the 
experimental values. 

3.5. Relation among < T >, T and r 

The primary observables in a microlensing experiment are 
the optical depth r, the rate T and the average duration of 
the events < T > . We saw that the duration of a particular 
microlensing event is T = Re/vi-± and its average, Eq. 
(^ij) , gives the time scale of the microlensing events 



<T>- 



<T^> _J Teii^i)dNey _ J Teiifi)dr 



< TO > 



/ £o{p)dN^ 



Jeo{fi)dr 



(43) 



Assuming £i(/i) = £o(/^) = constant, we have that the 
time scale is then equal to < T >= {J T dr)/r. Recalling 
Eqs.(0) and (|3C|), the relation among the three observables 
is just 



<r>=— — ^</i >, 



7i"r pq. 



lens 



(44) 



where we have set 



ratio 



M, 



1. If /iup is equal to ^max, the 



© 



PO.lo 



< /i^ >= 1 (see Eq. (EQ)), and we get the 



known relation < T >— 



4. Estimate of microlensing quantities for LMC 
self-lensing 

With the formalism developed in the previous sections, we 
are able to get the theoretical estimates for several impor- 
tant microlensing quantities. These estimates are obtained 
using the values for the parameters as described in Sect. 
P (only the velocity dispersion and the experimental effi- 
ciency are left unspecified, although we will also consider 
the influence on our results of these parameters). While 
the uncertainties of the scale lengths and the LMC disk 
inclination induce only minor modifications in our results, 
the larger uncertainties on the LMC mass are much more 
important. All results are summarized in Tables ^ and |^ 
for the models la and 16, and in Tables ^ and | for the 
models 2a and 2b. The LMC mass (without the halo) is in 
the range 3 — 6 x 10^ Mq. Four different IMFs were used 
as specified below. We notice also that using the above 
discussed source and lens distributions for the LMC we 
did not take into account that there are also gas and dust 
clouds which will obscure some regions of the LMC and 
thus due to that the values for the optical depth and espe- 
cially for the expected number of self-lensing events have 
to be considered as upper limits. 
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Galactic latitude 
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Fig. 1. The 27 MACHO fields that we selected. The filled 
stars denote the positions of the microlensing events re- 
ported by the MACHO collaboration, while the empty 
stars are that of the EROS collaboration. The centers of 
the LMC bar and disk are labeled by "b" and "d" , respec- 
tively. The blue line marks the LMC bar major axis. 

4.1. Optical depth 

We computed using Eq. the optical depth for several 
different source/lens geometries. For each of them, adopt- 
ing the parameters discussed in Sect. ^ we calculated the 
optical depth in the directions of 27 MACHO fields, see 
Fig. |l|. These fields (1, 2, 3, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 
17, 18, 19, 22, 23, 24, 47, 76, 77, 78, 79, 80, 81, 82 accord- 
ing to the MACHO numbering) have been selected taking 
into account the position of the microlensing events found 
towards the LMC by the MACHO collaboration during a 
5.7 years monitoring campaign (Alcock et al. 2000a). 

Following the strategy of Gyuk et al. (2000), the values 
of the optical depths, found for each source/lens geometry, 
have been averaged on the 27 MACHO fields that we se- 
lected. The value for each field was obtained by computing 
the optical depth per unit surface and then by integrating 
over the area of the field. 

• Model 1. The total optical depth, obtained using the 
"preferred" values for Mdisk + Mbar = 3 x 10^ Mq and 
an inclination angle of 30°, is r = 2.63 x lO^'* for the 
model la and r = 2.58 x 10"** for the model 16. The 
results, summarized in Table ^ are in good agreement 
with the corresponding values found by Gyuk et al. 
(2000) as expected, since we adopted the same values 
for the LMC parameterization. There is a slight depen- 
dence on the LMC disk inclination, as it appears also 
from Fig. ^, which shows the variation of the total opti- 
cal depth with the disk inclination. If we use the upper 
limit for the LMC disk and bar masses, the correspond- 
ing total optical depth is clearly larger: r = 5.25 x 10~^ 
for the model la and t = 5.15 x 10~® for the model 



3.4 



3.2 



- 3 



2.8 



2.6 



Model 1-a 



Model 1-b 



30 32 34 36 

LMC inclination angle 



38 



40 



Fig. 2. Total optical depth for LMC self-lensing as a func- 
tion of the LMC inclination angle for two different models 
assuming a total mass of 3 x 10^ Mq. 

16. By appropriately choosing the parameters in the 
ranges specified in Sect. 2 Gyuk et al. (2000) found an 
upper value of about r = 8 x 10"^ for the self-lensing 
optical depth. Clearly, this implies to adopt rather un- 
realistic values for the parameters. We refer to Gyuk et 
al. (2000) for a detailed discussion on the self-lensing 
optical depth range as well as a comparison with the 
theoretical results appeared meanwhile in the litera- 
ture. 

• Model 2. The total optical depth is r = 3.46 x 10~® for 
the model 2a and t = 4.05 x 10~® for the model 26 for 
the lower mass value. These values are quite similar to 
the ones found for model 1. The adopted parameters 
are described in Sect. |^ and the results are summa- 
rized in Table |^ for different source/lens geometries. 
As in the previous case, if we use the upper limit for 
the LMC disk and bar masses, the corresponding total 
optical depth gets larger and probably unrealistic for 
both models. Nonetheless, even in that case the ob- 
served value for the optical depth cannot be explained 
only by self-lensing. 

The five events found so far by the EROS collaboration 
are all located quite far from the LMC central regions, see 
Fig. |l|. Therefore, we expect the optical depth for self- 
lensing to be accordingly small. Indeed, using the model 
la, the averaged self-lensing optical depth over the event 
positions is 0.8 x 10~* for a LMC disk mass equal to 2.4 x 
10^ M0 or 1.6 X 10"^ for a LMC disk mass equal to 4.8 x 
10^ Mq (given the position of the EROS events the bar 
does not contribute) . Clearly, these values are smaller than 
the corresponding ones for the MACHO fields and thus the 
EROS events are less likely to be due to self-lensing. 

In Table ^ we give also the value for the optical depth 
due to MACHOs located in the LMC halo assuming a 
core radius a equal to 2 kpc (it was found that the optical 
depth is quite insensitive to the variation of the core radius 
(Gyuk et al. 2000)) and that MACHOs make up 100% of 
the dark LMC halo, which is clearly an extreme assump- 
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Table 1. Model 1: theoretical estimates of the optical depth, r, towards the LMC, averaged over 27 MACHO fields, 
for different source/lens geometries and parameters in the case of self-lensing. 



SouTce j Lens 


LAIC mass 


Model 


LMC 


Relative 


T 


Geometvy 






inclinatioTi 


W eight 








la 


30° 


0.40 


1.27 - 2.54 






la 


35° 


0.35 


1.38 - 2.76 


Disk/ disk 




la 


40° 


0.31 


1 52 - 3 04 


3-6 


lb 


30° 


0.69 


1.27 - 2.54 






lb 


35° 


0.64 


1.38 - 2.76 






lb 


40° 


0.60 


1.52 - 3.04 






la 


30° 


0.40 


83 - 1 66 






la 


35° 


0.35 


0.74 - 1.47 


Disk/bar 




la 


40° 


0.31 


66 - 1 32 


3-6 


lb 


30° 


0.69 


80 - 1 5Q 






lb 


35° 


0.64 


0.78 - 1.56 






lb 


40° 


0.60 


0.71 - 1.42 






la 


30° 


0.60 


1 20 - 2 40 






la 


35° 


0.65 


1.50 - 2.94 


Bar /disk 




la 


40° 


0.69 


1 79 - 3 59 

-i.it/ T J . t_i f_y 


3-6 


lb 


30° 


0.31 


2 31 - 4 62 






lb 


35° 


0.36 


2.76 - 5.53 






lb 


40° 


0.40 


3.29 - 6.57 






la 


30° 


0.60 


1 78 - 3 58 






la 


35° 


0.65 


1.91 - 3.84 


Bar /bar 




la 


40° 


0.69 


2.21 - 4.43 


3-6 


lb 


30° 


0.31 


1 38 - 2 76 






lb 


35° 


0.36 


1.50 - 3.00 






lb 


40° 


0.40 


1.65 - 3.30 






la 


30° 


1 


2.63 - 5.25 






la 


35° 


1 


2.99 - 5.88 


Total 




la 


40° 


1 


3.44 - 6.88 


3-6 






lb 


30° 


1 


2.58 - 5.15 






lb 


35° 


1 


2.92 - 5.84 






lb 


40° 


1 


3.33 - 6.66 



Table 2. Model 2: theoretical estimates of the optical depth, r, towards the LMC, averaged over 27 MACHO fields, 
for different source/lens geometries and parameters in the case of self-lensing. 



Source/ Lens 


Model 


LMC mass 


Relative 


r 


Geometry 




(xIO'^Mq) 


Weight 


(xlO^*^) 


Disk / disk 


2a 
2b 


3 - 6 
3 - 5 


0.50 
0.20 


2.14 - 4.61 
1.34 - 2.40 


Disk/bar 


2a 


3 - 6 


0.50 


1.79 - 3.59 


2b 


3 - 5 


0.20 


4.48 - 7.47 


Bar' / dink 


2a 
2b 


3 - 6 
3 - 5 


0.50 

o.cso 


1.31 - 2.63 
0.82 - i.:i6 


Bar /bar 


2a 


3 - 6 


0.50 


1.68 - 3.36 


2b 


3 - 5 


0.80 


2.79 - 6.99 


Total 


2a 


3 - 6 


1 


3.46 - 7.10 


2b 


3 - 5 


1 


4.05 - 8.65 



tion and has thus to be considered as an upper limit given 
also the fact that it is still not clear if the LMC has a 
halo at all. Moreover, this value is valid for a line of sight 
going through the center of the LMC and, therefore, also 
for that reason it has to be considered as an upper limit. 



4.2. Number of events 

By observing roughly 12 millions stars from 1992 to 
1998, the MACHO collaboration found 13-17 microlens- 
ing events, according to the data-analysis criteria used {A 
or B) (Alcock et al. 2000a). Instead, the ER0S2 collabora- 
tion found 4 events monitoring about 25.5 millions of stars 
from 1996 to 1999 (Milsztajn & Lasserre 2001). Another 
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Fig. 3. Bar chart of the number of observed microlensing 
events towards the LMC as a function of their duration. 
The full bars represent the seventeen MACHO events, 
while the gray bars represent the five EROS events. 

event was found by EROSl in the period 1990-1993 mon- 
itoring roughly 8 milion of stars (Aubourg et al. 1993). 
These results are summarized in the bar chart of Fig. 
from which one sees that the detection of events with short 
duration {i < 60 days) is favoured. Yet, the two MACHO 
events, LMC-9 and LMC-22, both with £ > 60 days, must 
be treated with caution, because they are a binary event 
(Alcock et al. 2000b) and a possible exotic event (Alcock 
et al. 2000a), respectively. 

We computed the number of self-lensing events, iVcv, 
again by averaging over the 27 MACHO field as done for 
the optical depth. We used also the MACHO parameters: 
iV* = 11.9 X 10^ stars, tots = 5.7years, e = 0.4 - 0.5. 
We recall that the number of events is N^v — N^,tohs^, 
where F is given by Eq. (^0[). For the mass distribution 
we used four different Salpeter-type IMFs as given in Eq. 
(p^), with, as an illustration, different values for a and 

i) IMF-1 with ^inin = 0.1 and a = 0.56 for 0.1 < ^ < 0.59 

and a = 2.2 for fi > 0.59; 
a) IMF-2 with /i,„i„ = 0.1 and a = 2.35; 
in) IMF-3 with /^mi„ = 0.01 and a = 0.8 for 0.01 < < 0.2 

and a = 2.35 for n > 0.2; 
iv) IMF-4 with ^„iin = 0.01 and a = 2.35. 

• Model 1. The results for the different source/lens ge- 
ometries, IMFs and parameters are shown in Table || 
for £ = 0.5, a = 30 km/sec. For other values of e and 
a it is sufficient to rescale the numbers in the table by 
the new values of the parameters, since the number of 
events is proportional to e and a. 
The IMF-1 is not very realistic since it is based upon 
counts of M class stars in the Galactic disk (Gould et 
al. 1997) and should thus be taken just as an illustra- 
tion. Moreover, as can be inferred from Table |, for the 
IMF-1 mass distribution the number of events is less 
than for the the more realistic IMF-2 and IMF-3 cases. 
The IMF-4 is also not very realistic, leading already 



for a rather small dispersion velocity to a conspicu- 
ous number of events. However, as we will see below, 
the resulting average event duration for the IMF-4 is 
definitely not in agreement with the observed value. 
For the IMF-2 and the IMF-3 we see that even by con- 
sidering the more extreme values for the various pa- 
rameters we get at most about 10 events, which is still 
less then the ones reported by the MACHO group al- 
though not very much. However, this implies to adopt 
rather unlikely values for the various LMC parameters. 
For the preferred and more realistic values, leading to 
an optical depth of t ~ 2 — 3 x 10~^, we expect about 
3 to 5 events due to self-lensing (Table ^ . 

• Model 2. For the model 2, on the ground of the results 
from the previous section, we concentrate only on the 
more realistic models for self-lensing by using IMF-2 
and IMF-3. The results for the different source/lens 
geometries are shown in Tabic || for model 2a and 26, 
where we fixed again the detection efficiency to 0.5 and 
the velocity dispersion to 30 km/sec. We see that also 
these models lead to values which are similar to the 
ones of model 1 and thus are also not able to explain 
all the observed MACHO events. For the model 2 we 
expect again rougly 3 events due to self-lensing when 
using the lower value for the LMC mass. Thus it seems 
that this conclusion is rather robust and does not de- 
pend much on the detailed form of the bar and the 
disk. 

Similarly, we can estimate the number of self-lensing 
events we expect for the ER0S2 data for a disk/disk 
source lens geometry. In this case, using the IMF-2 and 
an efficiency equal to about 0.12 (see Milsztajn & Lasserre 
2001), the estimated number of events, using the model la, 
is 0.16 — 0.32 for cr ~ 30 km/sec and for a LMC disk mass 
in between 2.4 - 4.8 x IO^Mq. For the IMF-3 the corre- 
sponding number of events is 0.19 — 0.37. If instead we use 
the model 2a, we expect for the above parameters a num- 
ber of events equal to 0.23 — 0.46 using the IMF-2, and 
0.27-0.53 using the IMF-3. Finally, we expect 0.14-0.24 
and 0.17 - 0.28 for the model 2b and for IMF-2 and IMF- 
3, respectively. Even if the small number of events found 
by ER0S2 does not allow to draw a definitive conclusion, 
it seems unlikely that the ER0S2 events are all due to 
self-lensing. 

We computed also the expected number of events due 
to MACHOs in the LMC halo. The lens number density 
is given by 

^=fl + ^)"^, (45) 
a/i \ a- / a/x 

where for dna/dfi we assume a delta function, and fix a 
equal to 2 kpc (Gyuk et al. 2000) (sec Table |). Using the 
MACHO collaboration parameters and ~ 0.5 we expect 
roughly 5 events under the assumption of a 100% MACHO 
LMC halo and assuming that the rate is uniformly equal 
to the one for the line of sight going through the center 
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Table 3. Model 1: theoretical estimates of the LMC self-lensing number of events, A'ev, averaged over 27 MACHO 
fields. We used e — 0.5 for the detection efficiency and a = 30kmsec~^ for the velocity dispersion. Different source/lens 
geometries and four Salpeter-type IMFs with different slopes and minimal lens mass as described in the text were 
used. 



Source/ Lens 


LMC mass 


Model 


LMC 


Relative 


iVcv 


iVov 




iVov 


Geometry 


(xlO'-'MfTi) 




inclination 


Weight 


(IMF - 1) 


(IMF - 2) 


(IMF - 3) 


(IMF - 4) 






la 


30° 


0.40 


0.60-1.20 


1.50-2.85 


1.62-3.23 


4.35-8.85 


Disk /disk 


3-6 


la 
lb 


35° 
30° 


0.35 
0.69 


0.60-1.35 
0.60-1.20 


1.50-3.00 
1.50-2.86 


1.72-3.44 
1.62-3.23 


4.65-9.45 
4.35-8.85 






lb 


35° 


0.64 


0.60-1.35 


1.50-3.00 


1.72-3.44 


4.65-9.45 






la 


30° 


0.40 


0.45-1.05 


1.20-2.25 


1.31-2.63 


3.60-7.20 


Disk/bar 


3-6 


la 
lb 


35° 
30° 


0.35 
0.69 


0.45-0.75 
0.15-0.30 


0.90-1.80 
0.30-0.60 


1.04-2.08 
0.33-0.66 


2.85-5.70 
0.90-1.80 






lb 


35° 


0.64 


0.15-0.15 


0.30-0.45 


0.30-0.60 


0.75-1.65 






la 


30° 


0.60 


0.90-1.80 


2.10-4.35 


2.50-5.00 


6.75-13.7 


Bar /disk 


3-6 


la 
lb 


35° 
30° 


0.65 
0.31 


1.35-2.55 
2.25-4.50 


3.00-6.00 
5.10-10.2 


3.47-6.94 
5.95-11.9 


9.45-19.1 
16.4-32.6 






lb 


35° 


0.36 


2.85-5.70 


6.60-13.1 


7.57-15.1 


20.7-41.4 






la 


30° 


0.60 


1.05-2.10 


2.40-4.80 


2.81-5.62 


7.65-15.5 


Bar /bar 




la 


35° 


0.65 


1.20-2.25 


2.55-5.35 


3.03-6.06 


8.25-16.7 


3-6 


lb 


30° 


0.31 


0.75-1.35 


1.50-3.15 


1.81-3.62 


4.95-9.90 






lb 


35° 


0.36 


0.75-1.50 


1.65-3.45 


1.97-3.93 


5.40-10.8 






la 


30° 


1 


1.65-3.30 


3.75-7.50 


4.36-8.72 


11.9-23.9 


Total 


3-6 


la 


35° 


1 


2.10-3.90 


4.35-9.00 


5.19-10.4 


14.1-28.5 


lb 


30° 


1 


1.50-2.85 


3.30-6.45 


3.75-7.50 


10.2-20.6 






lb 


35° 


1 


1.80-3.60 


4.05-8.10 


4.72-9.44 


12.9-25.9 



Table 4. Model 2: estimates of the total number of microlensing events towards LMC for the MACHO observation 
time computed in the case of self-lensing and for different values of the velocity dispersion, LMC mass and efficiency. 
We used e = 0.5 for the detection efficiency and a = 30 km see" ^ for the velocity dispersion. Two Salpeter-type IMFs 
were used (see text for the values of the parameters). 



Source/ Lens 


LMC 


mass 


Model 


Relative 


iVcv 




Geometry 


(xlO^ 


Mq) 




Weight 


(IMF 


-2) 


(IMF 


-3) 


Disk/disk 


3 - 


6 


2a 


0.50 


2.38 - 


4.77 


2.76 - 


5.53 


3 - 


5 


2b 


0.20 


1.49 - 


2.48 


1.73 - 


2.88 


Disk /bar 


3 - 


6 


2a 


0.50 


0.58 - 


1.17 


0.68 - 


1.35 


3 - 


5 


2b 


0.29 


1.46 - 


2.43 


1.69 - 


2.82 


Bar /disk 


3 - 


6 


2a 


0.50 


2.11 - 


4.23 


2.45 - 


4.91 


3 - 


5 


2b 


0.80 


1.32 - 


2.21 


1.53 - 


2.56 


Bar /bar 


3 - 


6 


2a 


0.50 


0.61 - 


1.22 


0.70 - 


1.41 


3 - 


5 


2b 


0.80 


1.52 - 


2.53 


1.76 - 


2.93 


Total 


3 - 


6 


2a 


1 


2.84 - 


5.70 


3.30 - 


6.60 


3 - 


5 


2b 


1 


2.86 - 


4.77 


3.32 - 


5.53 



of the LMC, this last assumption obviously leads to an 
upper limit. Clearly, since the halo mass fraction in form 
of MACHOs in the halo of our Galaxy is at most 20% 
(Alcock et al. 2000a), we do not expect this fraction to 
be higher in the LMC halo. Thus, realistically we expect 
at most about 1 event coming from this component. A 
similar conclusion holds for the ER0S2 data. 

4.3. Average event duration 

Using the Eq. ( pl| ) we get estimates of the average event 
duration < T > for different source/lens geometries and 
IMFs. The results, averaged on 27 MACHO field, are 



shown in Table || for model 1 and in Table for model 
2. The velocity dispersion is in both cases assumed to be 
30 km/sec. 

• Model 1. In Fig. |^ we report the average event duration 
obtained using the four considered IMFs and different 
values of the velocity dispersions. We recall that the 
IMF-1 was deduced by dwarf star counts in the Milky 
Way disk, so we must be cautious when applying it to 
a population of another galaxy like LMC. 
The results in Fig. ^ can not be directly compared 
with the observed MACHO event durations as sum- 
marized in Table ^ because among them there are 
two events, LMC-9 (binary event) and LMC-22 (ex- 
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Table 5. Model 1: theoretical estimates of the LMC self-lensing event durations, < T >, averaged over 27 MACHO 
fields. We used a — SOkmsec^^ for the velocity dispersion. Different source/lens geometries and four Salpeter-type 
IMFs were used (For the parameters see text). The durations are given in days. 



OUlLi Cfc./ J^fc./fcL» 


IVl Ucifc-t 


T Mn 

IjlVl o 


^ T '--^ 


<^ T ^ 


<^ T ^ 




GcoTTiGtvy 




iTLclindtioTi 


(IMF — 1) 


{IMF — 2) 


{IMF — 3) 


{IMF — 4) 




1 a 

_LcL 


ou 


u / .o 




42.9 


1 Q q 
1 y . y 


Disk /disk 


let 


OO 


Do.o 


Pin 7 
ou. ( 




90 


1 h 


ou 


U i .o 


4:y.o 


49 Q 


1 Q Q 
±y.y 




lb 


35° 


68.8 


50.7 


43.9 


20.3 




Td 


ou 


Oi/.V 


A'^ Q 
4:0. y 


oo.u 


1 7 fi 


Disk /bar 


1 Q 

Id, 


00 


DU. / 


/l/l 7 


OO. 1 


1 7 Q 


1 h 


ou 




Af\ 7 




1 8 7 




lb 


35° 


64.3 


47.3 


41.0 


19.0 




la 


ou 






47 Q 


22.2 


Bar /disk 


la 


00 




oy.z 


01. o 


ZO. ( 


lb 


OU 


OZr.O 


(SO fi 


^9 ^ 


24.3 




lb 


35° 


87.1 


64.1 


55.5 


25.7 




la 


30° 


54.3 


40.0 


34.6 


16.0 


Bar /bar 


la 


35° 


54.3 


40.0 


34.6 


16.0 


lb 


30° 


65.4 


48.1 


41.7 


19.3 




lb 


35° 


65.3 


48.1 


41.7 


19.3 




la 


30° 


64.1 


47.2 


40.9 


18.9 


Average 


la 


35° 


66.1 


48.7 


42.1 


19.5 


lb 


30° 


69.6 


51.2 


44.4 


20.6 




lb 


35° 


71.4 


52.6 


45.5 


21.1 



Table 6. EROS and MACHO (blended and unblended) 
average event durations. The quantity tst{A, B) is the av- 
erage event timescale, according to the two different crite- 
ria used in the efficiency, for events in the MACHO Monte 
Carlo calculations that have been detected with an un- 
blended fit timescale of t. 



Team 


Set 


iVcv 


< T > 


EROS 


t 


5 


33 days 


MACHO 


i 


17 


44 days 


MACHO 


ist{A) 


13 


47 days 


MACHO 


ist{B) 


17 


56 days 



otic event), which are not compatible with our method 
of analysis and there is one more event, LMC-5 (disk 
event) which is not a self-lensing event. 
So, we must exclude these three events from the 
MACHO blended average event durations, set t, before 
comparing with our results. The remaining 14 events 
give an average duration of 37 days. Moreover, as men- 
tioned above (see also Fig. ^) all relevant events have 
a duration in between 18 and 52 days. Comparing this 
with our results (Fig. |^) and taking also into account 
that not all events are possibly due to self-lensing, 
clearly leads already to restrictions on the IMF and 
in particular on the velocity dispersion values: the ve- 
locity dispersion is more likely to be in the range of 
30 — 35 km/sec, that suggests a predominance of lenses 
of old population (Hughes et al. 1991) and IMF-2 or 
IMF-3 are favored, since lMF-1 leads to too high du- 
rations. On the other hand an IMF-4 model would re- 



Table 7. Model 2: theoretical estimates of the LMC self- 
lensing durations of events, < T >, averaged over 27 
MACHO fields. We used cr = 30 km sec" ^ for the veloc- 
ity dispersion. Various source/lens geometries and three 
Salpeter-type IMFs were used (for the parameters see 
text). The durations are given in days. (We do not dis- 
tinguish between model 2a and 2b, since both give the 
same results.) 



Source/ Lens 
Geometry 


<T> 
{IMF - 2) 


< T > 
{IMF - 3) 


< T > 
{IMF - 4) 


Disk / disk 


52.6 


45.6 


21.1 


Disk/bar 


46.3 


40.1 


18.5 


Bar /disk 


33.8 


29.3 


13.6 


Bar /bar 


41.6 


36.0 


16.7 


Average 


43.6 


37.8 


17.5 



quire a low velocity dispersion (below 20 km/sec), oth- 
erwise the average duration gets substantially below 
the observed value. However, for very low mass objects 
one would rather expect a higher velocity dispersion. 

• Model 2. 

For this model we consider three IMFs: IMF-2, IMF-3 
and IMF-4. The results, reported in Table ^ and in Fig. 
H, show that the mean durations are slightly lower than 
that of model 1, nonetheless the previous conclusions 
hold. 

We arrive at the same conclusion also if we compare our 
results with the EROS mean event duration (see Table |6|) 
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Fig. 4. Model 1: variation of the total self-lensing event 
durations as a function of the velocity dispersion. The solid 
lines and the dashed lines are obtained by using a LMC 
disk inclination of 30° and 35°, respectively. Four different 
IMFs were used (for the parameters see text). Two models 
(lo and 16) are shown separately. 



under the rather unplausible hypothesis that they are all 
due to self-lensing. 

For the IMF-4 mass distribution we see from Fig. || 
that the duration of events decreases slowly by increasing 
the velocity dispersion of the lenses and in order to explain 
the observed mean duration of the events requires a low 
(a < 20 km/sec) velocity dispersion. For what we saw in 
the previous section on the IMF-4 model this would then 
also imply that all the detected events are due to self- 
lensing. However, such a low required dispersion velocity 
is rather unplausible considering also the small mass of 
the lenses. Thus the IMF-4 model is certainly not very 
realistic. 

Using the number density given by Eq. (p5|), we esti- 
mated the event duration also for lenses in the LMC halo, 
see Table We expect a duration in between 128^^ and 
181^^ days according if we take into account the motion 
of the source as in the former case or not. 
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Fig. 5. Model 2: variation of the total self-lensing event 
duration as a function of the velocity dispersion. These 
curves are obtained by using a LMC disk inclination of 
35°. Three different IMFs were used (for the parameters 
see text). 



4.4. Self-lensing mean mass 

With the mass moment method we are able to get an esti- 
mate of the lens mean mass < M > , as given by Eq. (^2|) . 
We notice that to get this way the mean mass we do clearly 
not need to assume a particular IMF distribution. For 
< > and < >, we used the Einstein radius cross- 
ing times measured by the MACHO collaboration for the 
17 events that they detected (Alcock et al. 2000a). Some 
of them were excluded in our LMC self-lensing analysis: 
the Galactic disk event LMC-5, the binary event LMC-9, 
and the exotic event LMC-22 (probably the source is a dis- 
tant active galaxy (Bennett 2001)). With this exclusion, 
the events LMC-1 and LMC- 15 have the shortest dura- 
tion, while the events LMC-7, LMC-13 and LMC-14 have 
the longest duration. Moreover, LMC-14 is a self-lensing 
event (Alcock et al. 2001b). In Section 4.2 we saw that 



the expected number of self-lensing events is at least 3. 
Assuming thus that among the observed events (not in- 
cluding the 3 excluded ones as explained above, one of 
which being also probably a self-lensing) just 3 are due to 
self-lensing, we can estimate a minimum and a maximum 
mean mass using the following 3 events: LMC-1, LMC-14, 
LMC-15 (mean duration: 28.5 days) for < A/,nin >, and 
LMC-7, LMC-13, LMC-14 (mean duration: 50.5 days) for 

< Mmax >■ 

In Fig. ^ we report the results for different values of 
the velocity dispersion obtained using models la, lb and 
2 (we do not distinguish between model 2a and 26, since 
both give the same results). Any other combination in the 
choice of the events (or in particular by assuming that 
more than 3 events are due to self-lensing) would lead to 
a higher value for < Mmin > and to a lower value for < 
A^max >■ Thus the values given here have to be considered 
as lower and upper limit, respectively. These results have 
a slight dependence from the LMC disk inclination, as it 
appears from Fig. ^ which shows also the variation of the 
mean mass for a disk inclination between 30° — 35°. 
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Fig. 6. Total mean minimal and maximal mass of the 
lenses in solar mass unity as a function of the velocity 
dispersion (km/sec) for models la, 16 and 2 (dotted line). 
For the models la and lb, the solid lines and the dashed 
lines are obtained by using a LMC disk inclination of 30° 
and 35°, respectively. 



We see that for our preferred values for the LMC mass 
and for <t ~ 30 — 35 kmsec"^, we get a mass in the range 
~ 0.1 — 0.6 Mq, which is consistent with the various as- 
sumptions we made. This result suggests that the lenses 
are low mass stars, which is certainly not surprising This 
can also be seen as a consistency check on the modelling 
of the LMC we adopted. 



5. Galactic and halo lensing populations 

In this section we turn to the case that some of the lenses 
are located in the Milky Way or in its halo. A Galactic 
model with three main components fits well all the dy- 
namical observations. These components are the disk, the 
spheroid and the halo. The halo, composed essentially of 
dark matter, contributes for most part of the total mass 
of the Galaxy and should give the higher contribution 
to the microlensing surveys. However, since at least one 
event (LMC-5) is due to a Galactic disk lens (Alcock et 
al. 2001c), we will examine two other important Galactic 
components: disk (thin and thick) and spheroid. In Tables 
and S we summarize the estimates of the microlensing 



quantities related to these components. The thin and thick 
disks are considered as separate components. 

There is no general agreement on the details of the 
structure of our Galaxy and in particular on the param- 
eters of its thick disk and spheroid. For the disk we vary 
the structural parameters in a wide range, while for the 
spheroid we refer to Giudice et al. (1994), where the mod- 
elling of this component is amply discussed. Moreover, we 
refer to Carrol & OstHe (1996), Roulet & MoUerach (1997) 
and Sparke & Gallagher (2000) for most of the parameters 
used to describe the Galactic components. 

5.1. Lenses in the Galactic Thick Disk 

We start our analysis on galactic components with the 
hypothesis that some lenses are associated with the old- 
est sub-population of the Galactic disk. We consider that 
these lenses, a fraction of which might be high velocity 
white dwarfs (Reid et al. 2001), are confined within the 
so-called thick disk. In this case, we used a lens number 
density distribution per unit mass dn/dfi given by 

dn ^ dnQ 
dfj, dfi 



exp 



X exp 



•\/i?o + -^oi ^os^ ^ ~ 2i?o£'oi cos b cos / 



i?o 



i?TD 



exp 



Rtd 
\Doi sin 5 1 



X (46) 



duQ 
dfi ' 



where Rq is the Earth distance from the center of the 
Galaxy, for which wc use the lAU recommended value of 
8.5 kpc. The common value for the thick disk scale length 
Rtd is 3.5 kpc (Bahcall 1986), even if its estimate varies in 
a wide range. Here, we will vary i?TD in between 2 and 5 
kpc. The thick disk has a larger scale height as compared 
to the thin disk. Its value is generally estimated to be in 
the range 1 < ztd < 1-5 kpc. For the local column den- 
sity we use Sj^ = 35 - 45 Mq pc'^ (Derue et al. 2001; 
Roulet & MoUerach 1997). All the estimates of microlens- 
ing quantities for Rtd = 3.5 kpc, and ztd = 1 (Bahcall 
1986), are reported in Tables || and ^. 

In Fig. ^ we plotted the variation of the optical depth 
as a function of -Rtd for different values of E™ and with 
values for ztd equal to 1 kpc (Fig. |^a) and to 1.5 kpc 
(Fig. 1^), respectively. From these figures one sees that 
the optical depth varies little with the scale length i?TD- 
Instead, there is a more clear dependence of r on the value 
of the scale height ztd, besides also the density. This is 
true also for the number and the mean duration of events, 
Fig. 0c-d and Fig. |^-f respectively, even if in these cases 
the parameters on which one has to vary are connected to 
the IMF and the dispersion velocity of the lenses. 

For i?TD — 3.5 kpc, ztd = 1, and S™ = 35 — 
45 Mq pc~^, we found that the optical depth towards the 
LMC is T = 3.09 - 3.97 x 10"^ (Table |), that is higher 
than the estimated white dwarf population optical depth 
T « 0.4 x 10-8 (Koopmans & Blandford 2001). 
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Table 8. Theoretical estimates of the optical depth for stars in the LMC due to different lensing populations: thin 
disk, thick disk, spheroid (stellar halo), MW halo and LMC halo. 



Lensing component 


Eo (Mqpc-") 


po (Mqpc-^) 


T 


Thin Disk 


(30 


-50) 




(0.92 - 1.54) X 10"** 


Thick Disk 


(35 


-45) 




(3.09 - 3.97) X 10"* 


Spheroid 






0.0015 X 10"^ 


4.29 X 10"** 


MWhalo 






0.0079 X 10"^ 


4.83 X lO"'^ 


LMC halo 






0.0223 X 10"^ 


8.36 X 10"* 



Table 9. Theoretical estimates of the number of events and event mean duration for LMC stars by different lensing 
populations: thin disk, thick disk, spheroid (stellar halo), MW halo and LMC halo. A ^-function was used to describe 
the mass distribution of the lenses in the halos. Instead, two IMFs with different slopes were used to describe the other 
Galactic components. Different values for the dispersion velocity, a, were used, e is the MACHO detection efficiency 
and /i is the lens mass in solar mass units. 



Component 


a 


dni) / d^ 


MACHO iVov 


EROS A^cv 


< T > 


Thin Disk 


(20 - 30)kmscc"^ 


IMF-1 
IMF-2 


(0.3 - 0.4) e 
(0.9-1.4) e 


0.04 - 0.05 
0.12 - 0.18 


(107-71) days 
(78 - 52) days 


Thick Disk 


(35 - 60)kmsec"^ 


IMF-1 
IMF-2 


(1.5-2.5) e 
(3.3 - 5.7) e 


0.2-0.3 
0.5-0.8 


(103- 60) days 
(76-44) days 


Spheroid 


(90 - 120) kmsec-i 


IMF-1 
IMF-2 


(2.8-3.8) e 
(6.5 - 8.7) e 


0.4-0.5 
0.9-1.2 


(69 - 52) days 
(51-38) days 


MW halo 


210kmsec"^ 


S 


107.7e^"^/^ 


14.6 M"'^' 


71 Y^days 


LMC halo 


50 km sec~^ 


5 


7.3eM"^/^ 


1.0^-^/^ 


(128- 181)^days 



Following the discussion in Sect. 3.3, for the mass dis- 
tribution we used two Salpeter type IMFs: IMF-1 with 
a = 0.56 for O.l < ^ < 0.59 and a = 2.2 for /i > 0.59; 
IMF-2 with a = 2.35. For MACHO experimental num- 
bers, in Fig. we report the variation of the number of 
events as a function of i?TD for different values of a and 
adopting e = 0.5 (Alcock et al. 2000a), ztd equal to 1 kpc 
(Fig. @c) and to 1.5 kpc (Fig. 0d). 

Here and in the next section, clearly the use of IMF-1 
is favoured, because it was obtained from Galactic disk 



observations, as we already mentioned in Sect. L2. In this 
case, joining thin and thick disk, we expect 1 — 3 events for 
Rtd =3.5 kpc, ztd = 1, and for the velocity dispersion 
in the range cttd = 35 — 60 km sec" ^ (Derue et al. 2001; 
Carrol & Ostlie 1996). 

At last, we estimated the mean duration of the events 
as a function of i?TD for different values of a and adopting 
ztd equal to 1 kpc (Fig. 0e) and to 1.5 kpc (Fig. Qf). 
It seems that a disk event was really observed, LMC-5 
(Alcock et al. 2001c). Its duration, 37.8 days, is not in good 
agreement with our predictions (Table |^), however one 
can not draw any conclusion on just one event. Moreover, 
we assumed a minimum mass in the IMF-1 of 0.1 Mq, 
whereas from the observation it turns out that the lens 
mass could well be below that limit, explaining thus its 
shorter duration. 

5.2. Lenses in the Galactic Thin Disk 

We consider now the possibility to have lenses in the 
thin disk. This component, contrary to the thick disk 
that is composed of low metallicity stars, is made of 



relatively young metal-rich stars with a small velocity 
dispersion trtd — 20 — 30 km sec" ^ (Roulet & MoUerach 
1997). We modelled it as in the previous case, see Eq. 
(^6|), but with a smaller scale height ztd = 0.1 — 0.35 kpc 
(Roulet & MoUerach 1997) and with a local column den- 
sity = 30 - 50Mqpc"2 (Roulet & MoUerach 1997; 
Derue et al 2001). 

For the thick disk we saw that the microlensing quanti- 
ties do not vary much as a function of the scale length, but 
there were variations with the scale height. For the thin 
disk the scale length is know with only a small range of un- 
certainty and, therefore, the modifications on our theoret- 
ical estimates will be negligible. Here we use Ztd = 0.3 kpc 
(Bahcall 1986), and the optical depth turns out to be 
T = 0.9 - 1.5 X 10"S, see Table | 

Using again two different IMFs as in the previous case, 
we estimated the number of events and their mean dura- 
tion. The results are reported in Table We see that the 
expected average event duration for both the thick and 
thin disk components is above 50 days when considering 
the more realistic IMF-1, so that it is rather unlikely that 
more than 1-2 events among the observed ones can be 
due to these populations, unless there are a lot of brown 
dwarfs, but these would have been discovered in the same 
way as in the case of the LMC-5 event. 

5.3. Lenses in the Galactic Spheroid 

Here we examine microlensing events by the stellar halo, 
i.e. the so-called Galactic spheroid. Looking at high lati- 
tudes far above the disk, there is a population of globu- 
lar clusters, old metal-poor stars or high-velocity stars. 
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Fig. 7. Galactic thick disk microlensing quantities as a function of the scale length, Rtb , for different local column 
densities, IMFs, dispersion velocities, and for the scale height, ztd, equal to 1 kpc (a, c, e) and 1.5 kpc (b, d, f). 
MACHO collaboration experimental values were used. Here we take the detection efficiency e equal to 0.5. (a, b) 
Optical depth: the solid hne was estimated for — 45Mqpc~^, the dashed line for E™ — 40Mqpc^^, and the 
dotted line for S™ = 35Mqpc~^. (c, d) Number of events, (e, f) Mean duration of events. 



supported by a large velocity dispersion, Ogp ~ (90 — 
120)kmsec-i (Carrol & Ostlie 1996; Roulet & MoUerach 
1997). The spheroid, approximately spherical, provides an 
important contribution to the size of the central part of 
the Milky Way, while its radial extension is smaller than 
that of the halo. Dynamical measurements and infrared 
surveys predict for the spheroid a mass roughly equal to 
5 — 7 X 10^°Mq, larger than that estimated on the basis of 
star counts. This suggests that the spheroid must be heavy 



and mostly dark (Caldwell & Ostriker 1981; Bahcall et al. 
1983). 

Giudice et al. (1994) have already discussed in some 
detail the possibility to have lenses in the spheroid. Using 
five different models based on the works of Caldwell & 
Ostriker (1981), Rohlfs and Kreitschmann (1988), Bahcall 
et al. (1983), they estimated the main microlensing quan- 
tities: an optical depth in the range (3.1 — 4.9) x 10~^, 
an event rate in the range (0.01 - 0.16)/^ an 



Ph. Jetzer, L. Mancini, G. Scarpetta: Microlensing towards the LMC 



17 



average event duration in the range (66 — ^fji days, 
where as usual, is the lens mass in solar mass units. 

Here we adopt the same simplified density profile, 
which fits well the Ostriker and Caldwell model as that 
in Roulet & MoUerach (1997), given by 



dn 



sp 



dfi 



(47) 



^Rq 



Dli-2RoDoi cosi9) 



1/4 



/Re 



where Rc — 0.17 kpc is the core radius (Roulet & 
MoUerach 1997) and — 82° is the angle between the 
line of sight of the LMC and the direction of the Galactic 
center. In this case, using a local density equal to = 
1.5 X 10"'^ Mq pc"'^, our estimate for the optical depth is 
T = 4.29 X 10~*, (see Table ^) which is in good agree- 
ment with the results of Giudice et al. (1994) and lies well 
within the range they found. 

Also in this case we used two different IMFs to es- 
timate the number of events and their mean duration. 
The results are reported in Table |^ from which we see 
that the average duration is somewhat higher than the 
observed value, unless the dispersion velocity is closer to 
120 km sec~^ (or if the IMF goes well into the brown dwarf 
region). In this latter case some (2 to 6 events) among 
the MACHO longer duration events, close to about 40-50 
days, could be due to lenses in the spheroid. This result 
is still in a good agreement with Giudice et al. (1994). 
Similarly, we expect at most 1 to 2 EROS events to be- 
long to the spheroid. Clearly, given the calculated mean 
event duration and number of expected events we see that 
it is well plausible that some of the observed events are 
due to lenses located in the spheroid, however, certainly 
not all of them. 



5.4. Lenses in the Galactic Halo 

At last we estimated the microlensing quantities in the 
hypothesis that the lenses are MACHOs located in the 
Galactic halo (Jetzer 1996). The MACHO number density 
distribution that we used is now 
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where d — 82° is again the angle between the line of sight 
and the direction of the Galactic center and R^ is the core 
radius of the dark matter, that is also uncertain. Bahcall 
et al. (1983) estimated Rc equal to 2 kpc, while Caldwell 
& Ostriker (1981) suggested 7.8 kpc. Alcock et al. (2000a) 
used instead 5 kpc. Here, following Primack et al. (1988), 
we use a core radius equal to 5.6 kpc. For dno/ dfi, we used 
a delta function, Eq. (|l|), with = 7.9 x 10"^ Mq pc"^. 

Fixing the usual halo value for the velocity dispersion, 
an = 210kmsec~^, we estimated the optical depth, the 



number of events and the mean time duration. The re- 
sults are tabulated in Tables || and ||. For fi = 0.5 Mq 
(as an illustration) we expect for a halo made entirely of 
MACHOs 61 events using MACHO efficiency e = 0.4, and 
76 events using MACHO efficiency e — 0.5, respectively. 
Similarly for ER0S2, where we use an efficiency equal to 
e = 0.12, we expect 19 events. The expected mean event 
duration is 50 days for fi — 0.5. 

For the optical depth we get the well-known value of 
~ 5 X 10"'', while the estimated optical depth for a halo 
white dwarf population is t = 0.13 x 10^^ (Koopmans & 
Blandford, 2001). 

6. Mean Mass 

From the previous discussion it is clear that almost cer- 
tainly not all the observed events (in particular among 
the ones detected by the MACHO collaboration) originate 
from the same lens population. Accordingly, it is difficult 
to estimate a mean value for the mass, since we do not 
know in general to which population each event belongs. 
In the following we will consider different possibilities of 
how to divide the events among the various populations 
and so get estimates of the masses using the mass moment 
method. Clearly this has to be considered as an illustra- 
tion, nonetheless it allows to get some insight and also 
to get some upper and lower limits which can already be 
useful. 

The five EROS events, due to their positions (see 
Fig. |l|) and durations, can hardly contain self-lensing and 
Galactic disk lenses. So most of them if not all (with per- 
haps the possibility of having 1 or 2 events due to the 
spheroid) are MACHOs in the Galactic halo, in which case 
we find the following mean mass: < M >= 0.26 Mq, as- 
suming a spherical halo model. If the 1 or 2 events with the 
longer duration are due to lenses in the galactic spheroid 
then of course we would get a smaller average mass for 
the remaining 3 to 4 halo events. 

The MACHO events contain most probably some 
lenses in the LMC and at least one event due to a Galactic 
component. We exclude as before from our analysis the 
MACHO events LMC-9 and LMC-22, whereas LMC-14, is 
possibly a self-lensing event. We will consider, as an illus- 
tration, four different models by dividing up the observed 
events among the various populations, we shall label them 
with: M-1, M-2, M-3 and M-4 (see below). 

We already saw in section 4.4 that we expect at least 3 
self-lensing events. For these three events, fixing the LMC 
dispersion velocity to 30 km/sec and clearly assuming that 
one of them is the LMC-14 event, we estimated a mini- 
mum (models M-1 and M-3 in which LMC-1, LMC-14 and 
LMC-15 are the self-lensing events) and a maximum (mod- 
els M-2 and M-4 in which LMC-7, LMC-13 and LMC-14 
are the self-lensing events) mean mass, see Table ^ and 
Fig. I 

If only the LMC-5 event is a disk event (we assume 
it to be in the thick disk), we found that its mean mass 
varies from 0.08 Mq for cttd = 35 km/sec (model T-1) to 
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Table 10. Theoretical estimate of the lens mean mass 
for microlensing towards LMC assuming the lens to be in 
different components. 
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0.11 Mq 


Self — lensing 


M-2 
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3 
3 


0.43 Mq 

0.11 Mq 




M-4 


3 


0.43 Mq 




M-1 


11 


0.37 Mq 


Halo 


M-2 

M-3 


11 

5 


0.25 Mq 
0.23 Mq 




M-4 


7 


0.17 Mq 




M-1 







Spheroid 


M-2 
M-3 



6 


0.30 Mq 




M-4 


4 


0.27 Mq 


Thick disk 


T-1 
T-2 


1 
1 


0.08 Mq 

0.24 Mq 



0. 24 .q for CTTD = 60 km/sec, (model T-2), see Tabic ^ 
(assuming instead a velocity of 20 km/sec as mentioned 
by the MACHO collaboration (Alcock et al. 2001c) we 
find with the mass moment method a mass of 0.03 Mq, 
which compares quite well with the value inferred by the 
MACHO team). 

For the remaining MACHO events we considered two 
situations. 

1. They are all due to lenses in the Galactic halo. The 
corresponding mean mass is: < M >= 0.37 Mq us- 
ing the events LMC-4, LMC-6, LMC-7, LMC-8, LMC- 
13, LMC-18, LMC-20, LMC-21, LMC-23, LMC-25, 
LMC-27 (model M-1) (mean duration: 39.4 days); 
< M >= 0.25 Mq using the events LMC-1, LMC-4, 
LMC-6, LMC-8, LMC-15, LMC-18, LMC-20, LMC-21, 
LMC-23, LMC-25, LMC-27 (model M-2) (mean dura- 
tion: 33.4 days). 

2. Otherwise, they are due in part to lenses in the 
spheroid and in part to lenses in the halo. We divided 
the events between these two component preferring for 
the spheroid the events with longer durations. We con- 
sidered again two situations by assuming that in one 
case 6 events are due to the spheroid and in the sec- 
ond case only 4. In the model M-3 we estimated the 
mean mass for the halo by using the events LMC-4. 
LMC-8, LMC-18, LMC-20 and LMC-27 (mean dura- 
tion: 30.9 days), while for the spheroid we used the 
events LMC-6, LMC-7, LMC-13, LMC-21, LMC-23 
and LMC-25 (mean duration: 46.5 days). In the model 
M-4 the mean mass for the halo was obtained using 
the events LMC-1, LMC-4, LMC-8, LMC-15, LMC- 
18, LMC-20 and LMC-27 (mean duration: 27.2 days), 
while for the spheroid we used LMC-6, LMC-21, LMC- 
23 and LMC-25 (mean duration: 40.1 days). AU the 
results are shown in Table ^ 

We see that for the MACHOs in the Galactic halo we 
find a mean mass in the range (0.17 — 0.37) Mq, which is 



Table 11. Fraction of the local dark mass density in form 
of MACHOs in the Galactic halo, which we get using 
MACHO and EROS data. 
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EROS 2 
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in agreement with the previous value of 0.26 Mq obtained 
from the EROS data. It has to be noticed that these values 
are for a standard spherical halo, thus if we consider halo 
models with a flattening or with anisotropics in the veloc- 
ity space the above values can also vary substantially and 
in particular decrease (De Paolis et al. 1996; Grenacher et 
al. 1999). 

For the LMC halo we expect a lens mean mass in the 
range (0.02 - 0.16)Mq. 

7. MACHO mass fraction in the Galactic halo 

Another important quantity to be determined is the frac- 
tion / of the local dark mass density detected in form 
of MACHOs in the Galactic halo, which is given by (De 
Riijula et al. 1991): 

f EE < // > ~ 127pc3 < //i > . (49) 
Po 

Using the values given by the MACHO collaboration for 
their 5.7 years data, we estimated / for each of the four 
models considered before and also for the ER0S2 data 
assuming that all their events are due to MACHOs in the 
halo. The results, obtained by assuming again a standard 
spherical halo model as in Sec. are reported in Table 
0- 

We see that for the models we considered the frac- 
tion of halo dark matter in form of MACHOs varies be- 
tween 5% and 13%. Moreover, when neglecting contribu- 
tions from the spheroid (models M-1 and M-2) the values 
we get from MACHO and EROS collaboration are prac- 
tically similar. Clearly, if some events are indeed due to 
the spheroid (as in models M-3 and M-4) we should then 
take this component also into account for the EROS data. 
In which case the fraction we find from the EROS exper- 
iment gets accordingly smaller and is then again in good 
agreement with the MACHO value. 

8. Asymmetry of the observed events in the LMC 

A striking feature which comes out by examining the po- 
sitions of the events found by the MACHO and EROS 
collaborations is that there is a clear asymmetry in the 
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spatial distribution of the events with respect to the bar 
major axis. The events are concentrated along the exten- 
sion of the bar and in the south-west side of LMC (Fig. 
|l|). In particular, from MACHO data we have 12 events 
located in the south-west side and 5 events located in the 
north-east side. The ratio is just 12/5 = 2.4. To study 
the nature of this asymmetry, we divided the 27 MACHO 
fields that we selected in two groups. In the first group 
there are the fields located above the bar major axis and in 
the second those below. Four fields (11, 12, 47, 79) due to 
their location were split among the two groups. We added 
the products obtained by multiplying the number of stars, 
iV^, by the observing time, Tobs, that is clearly propor- 
tional to the number of observations pro field reported by 
the MACHO collaboration (Alcock et al. 2000a), weighted 
by the corresponding event rate for self-lensing F for each 
of the fields in one group. The ratio of the two sums turns 
to be close to 1, more precisely 

( Tohs,i Fi ) 

-j^^ ^ 0.96. (50) 

\ / Tie 

This result is obtained using the model la and IMF-2, 
although we would get a similar value also by using the 
other models. If the events are due to self-lensing we would 
expect the ratio of the events to be close to 1 rather than 
being 2.4. Clearly, the same holds if the events are due to 
MACHOs in the halo or in one of the considered galac- 
tic populations. A possible explanation might be that the 
lenses are located in the halo of the LMC which would 
not be spherical and perhaps more elongated towards our 
galaxy, or an alternative is that the MACHOs in our halo 
are not uniformly distributed. 

9. Discussion 

As mentioned at the beginning the issue of the location 
and nature of the objects which act as lenses in the ob- 
served microlensing events is still an open problem, which 
can possibly be solved once more events will be avail- 
able. The number of events found by the experimental 
collaborations are too few in comparison with that pre- 
dicted for a halo composed entirely by MACHOs. In the 
last years several possible explanations of the experimen- 
tal data have been proposed, but they are all not definitive 
and not completely satisfactory. One possibility might be 
that there are more MACHOs in the halo, but that, for 
instance, they are associated with gas cloud (Bozza et al. 
2002), which would then produce non-achromatic events, 
that due to the present selection criteria have not been 
considered. On the other hand, one should also consider 
the possibility that some events might not be due to mi- 
crolensing at all. This obviously underlines the fact that 
the present results have to be taken with care and that 
more observations are needed in order resolve this issue. 

From the presently few available events and our above 
discussion it emerges clearly that, especially for the 
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MACHO events, the lenses are due to different popula- 
tions. Some are certainly due to LMC self-lensing, but 
this can hardly be the case for all the observed events 
especially for the few EROS events. 

For the preferred LMC model with a dispersion veloc- 
ity of about 30 kmsec"^ we expect some 3 events among 
the MACHO ones (not including the possible self-lensing 
binary event LMC-9) to be due to self-lensing. Moreover, 
with the mass moment method we find an average mass in 
the range (0.1 — 0.5) Mq for the self-lensing events, which 
is consistent with the expectation that the lenses are low 
mass stars. The contribution of a LMC halo is, even if it 
exists, also very minor of at most 1 event unless we are 
in the rather strange situation where the LMC halo is, 
contrary to our own, made almost entirely of MACHOs. 
From the galactic component, thin and thick disk, we ex- 
pect roughly 1 — 3 events on the MACHO data, but no 
contribution on the EROS data. This result seems to be 
in agreement with the fact that LMC-5 is a disk event. 
The inferred mass is small, though compatible with a low 
mass star, but clearly with only one event at disposal one 
has to take this value just as indicative. 

As a result of our analysis we find that a plausible so- 
lution is that among the MACHO data some 3 — 4 events 
are due to self-lensing, 1 — 2 to the thick disk and for the 
remaining ones about half are due to the spheroid and 
the others to the halo. If so the optical depth gets con- 
tributions from each of these components, namely: about 
2.3 X 10~* from self-lensing, some (2 — 4) x 10~^ from the 
halo, 4x 10~* from the spheroid and 4x 10~^ from thin and 
thick disk. This way we get a total optical depth of about 
(1.2—1.5) X 10~^ which is in good agreement with inferred 
value by the MACHO team of r = 1.21°;^ x 10"^. Since 
the EROS events are less and given their spatial position, 
we do not expect that they get much contribution from 
self-lensing and the thick disk so that it is clearly not pos- 
sible to draw more conclusions from them. On the other 
hand assuming that the EROS events are due to lenses 
in the halo, with or without some contribution from the 
spheroid, leads to a halo mass fraction which is in rea- 
sonable agreement with the corresponding MACHO value 
once the self-lensing and the disk events are subtracted. 
This way the MACHO and EROS results nicely fit to- 
gether. Moreover, the above mentioned asymmetry, if not 
just due to a statistical fluctuation, is also not compatible 
with only self-lensing events. Given our results it is also 
clear that one has to take with care values on the lens mass 
based on the assumption that all lenses belong to just one 
population. Clearly, once more data will be available, by 
using also the methods outlined in this paper, it will be 
possible to draw more firm conclusions. 
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Appendix A: Coordinate transformations 

To quantify pd or pb, we introduce the coordinate sys- 
tem {xq, Uo, zo} which has the origin in the center of the 
LMC at {l,b,D} = {lQ,bo, Dq}, and has Zq-slxis toward 
the observer, Xo-a.xis anti parallel to the galactic longi- 
tude axis, and j/o-axis parallel to the galactic latitude axis. 
Dq ~ 50 kpc is the distance between the center of the LMC 
and the observer, while D is the generic observer-lens or 
observer-source distance. {^0,^0} are the galactic coordi- 
nates of the center of the LMC. The coordinate transfor- 
mations are given by 

xo = —D cosbsm{l — Iq) 

yn = Dsinbcosbo — Dcosbsmbocos{l ~ lo) (A.l) 
zq = Df) — D cos b cos 60 cos {I — lo) ~ D sin b sin b^. 

The coordinate system of the LMC disk {x,y,z} is the 
same orthogonal system as {xo,yo, zq}, except that it is 
rotated around zg-axis by the position angle (j) counter- 
clockwise and around the new cc-axis by the inclination an- 
gle i clockwise. The coordinate transformations are given 

by 

X — Xo COS <p + yo sin 

y — —xosinipcosi + yQCOS(j)COsi — zosini (A. 2) 

z = —Xq sin (f) sin i + yo cos (j) sin i + zo cos «. 
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